A Spin- Orbital Singlet and Quantum Critical Point on the Diamond Lattice: FeSc2S4 
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We present a theory of spin and orbital physics in the A-site spinel compound FeSc2S4, which experimen- 
tally exhibits a broad "spin-orbital liquid" regime. A spin-orbital Hamiltonian is derived from a combination 
of microscopic consideration and symmetry analysis. We demonstrate a keen competition between spin-orbit 
interactions, which favor formation of a local "Spin-Orbital Singlet" (SOS), and exchange, which favors mag- 
netic and orbital ordering. Separating the SOS from the ordered state is a quantum critical point (QCP). We 
argue that FeSc2S4 is close to this QCP on the SOS side. The full phase diagram of the model includes a 
commensurate-incommensurate transition within the ordered phase. A variety of comparison to and suggestion 
for experiments are discussed. 
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The search for quantum spin liquids - materials in which 
local moments are well formed but continue to fluctuate quan- 
tum mechanically down to zero temperature - is a funda- 
mental challenge in condensed matter physics. In transition 
metal compounds with relatively high (e.g. cubic) symme- 
try, a richer possibility has also been contemplated, in which 
not only spin but also orbital states of localized electrons fluc- 
tuate. Such a "Spin Orbital Liquid" (SOL) was proposed in 
LaTiOs LU- Probably the best candidate for a SOL is FeSc2S4 
a spinel compound (with the general structure AB2X4), in 
which only the A sites form a magnetic ally/orbitally active 
diamond sublattice. In recent years, a variety of such A-site 
spinels, e.g. C0AI2O4 and MnSc2S4, were also found to be 
frustrated 12] [3] |4] |5] |6l, forming a "spiral spin liquid" |7| 
at certain temperature range. FeSc2S4 stands out markedly 
amongst this class of compounds in exhibiting a much broader 
liquid regime, extending down to the lowest measured temper- 
atures. 
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FIG. 1: (color online). The schematic phase diagram versus T, 
temperature and x, the ratio of exchange to spin orbit interaction 
(see text). Here the labels for the phases are: "SOS", spin-orbital- 
singlet; "CAF-l-OO", commensurate antifeiTomagnet with orbital or- 
der; "IC", incommensurate spin and orbital order. "QC" stands for 
the quantum critical regime. 



In the letter we describe a theoretical study of the spin and 
orbital physics of FeSc2S4 E 12 [TOl [H] [El . Here, the mag- 
netic Fe^+ ion is in a 'dd^ configuration, with a local S — 2 
moment as well as a two-fold orbital degeneracy (associated 
with one hole in the Cg doublet), i.e. Fe^+ is Jahn-Teller ac- 
tive. However, no indication of long range spin or orbital or- 
der is detected down to 50mK, which makes the frustration 
parameter / > 1000 (/ = Qcw /Tn and Qcw = -45. IK is 
the Curie- Weiss temperature), one of the largest values ever 
observed. We argue that the unusual behavior of FeSc2S4 
arises from proximity to a Quantum Critical Point (QCP) be- 
tween a spin-orbital singlet phase induced by atomic spin- 
orbit coupling and a magnetically and orbitally ordered phase 
favored by exchange (see Fig.[T]l. 

The spin-orbital Hamiltonian in FeSc2Si — The naively 
large (10-fold) spin and orbital degeneracy of the isolated 
Fe^+ ions must be split by additional effects. One natu- 
ral mechanism is exchange. To study this problem, we will 
first derive a "Kugel-Khomskii" type spin-orbital exchange 
Hamiltonian |13 |, in which the spin and orbital state on each 
Fe^+ are described by spin {S — 2) and pseudo-spin 
Ti (r — 1/2). The study of a structurally identical A- 
site spinel compound MnSc2S4 (no orbital degeneracy) in- 
dicates El [H [15] fliat a minimal superexchange Hamilto- 
nian should include both the nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) interactions. Unfortunately, direct 
derivation of the superexchange Hamiltonian from the "par- 
ent" Hubbard model is not very practical due to the com- 
plicated exchange paths A-X-B-X-A fTSl which involve five 
ions for both NNs and NNNs. However, the Fd3m space 
group symmetry of FeSc2S4 rather strongly constrains the 
exchange Hamiltonian. One can show that the most general 
spin-orbital exchange Hamiltonian allowed by symmetry for 
NNs and NNNs has the following simple form (neglecting the 
spin anisotropic terms), 

'Hex = \ y^\JijSi ■ Sj + ^KjjTj ■ Tj + kjjT^T^ 
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where Jij = Ji or J2 when ij are NN and NNN sites, respec- 
tively (and similarly for Kij, Kij, Lij, Lij). 

To further constrain the couplings, we must treat the micro- 
scopic physics. Integrating out the intermediate states on B 
and X sites, we consider processes involving electron motion 
between two Fe^+ ions. In general, superexchange may oc- 
cur between electrons from either the or eg levels. Only 
processes involving the latter, however, distinguish the orbital 
pseudo-spin. Now note that in Eq. ([T]), there are no terms lin- 
ear in the pseudo-spin operators t. Therefore, microscopic 
processes which might individually contribute to such terms 
for a single pair of Fe^+ ions must cancel when summed to- 
gether. Orbital contributions from superexchange processes 
between Cg and <2g are precisely of this form, and can there- 
fore be neglected. Thus to obtain the orbital part of the ex- 
change Hamitonian, we only need to focus on terms involving 
motion of two eg electrons. Again, symmetry restricts the 
form of Hamiltonian describing the effective transfer between 
NNs and NNNs: 

ij ma 

where d\^^^ and dima are the creation and annihilation oper- 
ators of a hole at mth eg level with spin cr, and = ti.t2 
when i^j are NNs and NNNs, respectively. Together with the 
on-site Hubbard-U terms, one can derive the pseudo-spin part 
of the exchange Hamiltonian by the standard perturbation the- 
ory. Combined with the contribution to the spin part from the 
t2g-t2g and t2g-eg exchange, we obtain a simplified version of 
Eq. ([T]l with 

Kij — Lij — 0, Lij — IKij, (3) 

with K12 = t\ 2/ (4?7) > 0. The simplified Hamiltonian con- 
tains 4 dimensionful parameters (Ji, J2, Ki, K2)- Moreover, 
from the above analysis, we expect crudely Ji 2 ^ K12 as 
the if 's only come from the Cg-eg exchange. 

The exchange Hamiltonian Ti^x has spin-rotational symme- 
try. A second mechanism to split the large ionic degeneracy 
is spin-orbit coupling. For an isolated Fe^+ ion, the physical 
spin and pseudo-spin can couple via the symmetry allowed 
term \TT\, 

K = -^{Vsrms^r - {sfr]+rn3{s!r ~ sm . (4) 

Note that, for an isolated single Fe^+ ion, Hq results in a non- 
degenerate ground state with a gap A to the first excited triplet. 
The ground state is a highly entangled state of spin and orbital 
degrees of freedom: a "spin orbital singlet". It is remarkable 
that a 3d ion with large spin S — 2 can form such a highly 
entangled quantum ground state. 

Such spin orbital singlet formation competes with ex- 
change, so it is helpful to have a microscopic estimate of A. 
This interaction arises at second order in the LS spin-orbit in- 
teraction Ao(L • S) due to the nonvanishing matrix elements 



of L between eg and t2g orbitals. Standard second order per- 
turbation theory gives [17J A = 6Ao^/Ate > 0, with Ate 
the crystal field splitting between eg and t2g levels. Taking 
the atomic spin-orbit coupling constant |Ao| ~ lOOcm"^ and 
Ate ~ 2500cm-i |18,i9J yields A w 36K. It is noteworthy 
that this is comparable to Qcw- If we assume Ki = K2 = 0, 
by the high temperature mean field theory (including both 
J^i'^h and Hex) one finds Ocw = -^^^(4Ji + 12 J2). 
Thus Qcw is a measure of the strength of exchange, and we 
conclude that exchange and spin-orbit couphng are competi- 
tive in FeSc2S4. 

From the above analysis, our complete Hamiltonian for 
FeSc2S4 is, in the first approximation 

^^E^O + ^ex. (5) 

i 

Minimal model for FeSc2Si — We begin the analysis by 
considering a simplified exchange Hamiltonian by appeal- 
ing to the neutron scattering measurements [gj, which indi- 
cate the low energy magnetic excitations are locaUzed near 
k = 27r(l, 0, 0). This is precisely the ordering wavevector as- 
sociated with a simple collinear Neel state on an FCC sublat- 
tice, and suggests the dominance of second neighbor exchange 
J2. Therefore we begin by studying the "minimal model" with 
Ji = Ki = K2 = and J2 > antiferromagnetic. The 
two FCC sublattices of the diamond lattice decouple in this 
case, and J2 can be viewed as a nearest-neighbor antiferro- 
magnetic exchange within either of these sublattices. (Patho- 
logical effects of this decoupling can be removed by including 
very small Ji). 

Here there is a single dimensionless parameter x = J2/A. 
For x 3> 1, the exchange is dominant, and since the S — 2 
spins are rather classical, we expect them to order magneti- 
cally. The ground states of the NN FCC sublattice are well- 
known. In real space, they consist of simple collinear antifer- 
romagnetic Neel states within each {100} plane, with an arbi- 
trary choice of axis for each such plane. In momentum space, 
this allows for spiral states with wavevector k = 27r(l, 5, 0) 
(and symmetry-related wavevectors) with arbitrary 5. 

The A term splits this degeneracy. For an arbitary mag- 
netically ordered state, in which (S^) ^ 0, the spin orbital 
Hamiltonian Eq. Q induces an "orbital field" that induces an 
orbital moment on each site. The magnitude of this orbital 
field is maximal when the spin is along one of the axial direc- 
tions [100]. This selects commensurate states with (5 = 0, 1/2. 
Within the minimal model, the remaining degeneracy is lifted 
by the weak effects of quantum fluctuations ll20l . which favors 
collinear spin states. Note that this selects 5 = Q, which corre- 
sponds to the experimentally observed low energy excitations 
in FeSc2S4. 

In contrast to the commensurate ordered phase (with 
collinear orbital order) found for a; ^ 1, for 2; ^ 1, the 
ground state is a spin orbital singlet, with a gap to all excited 
states. This is a generalization of a "quantum paramagnet" 
discussed intensively in spin-only models. The gap decreases 
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steadily upon reducing x, and is expected to close at a Quan- 
tum Critical Point (QCP). 

Indeed, this expectation is confirmed by a simple mean 
field theory (MFT). This consists of decoupling the exchange 
term as usual into an effective Zeeman field which is self- 
consistently determined for each site. Note that this procedure 
involves no approximation for TYq. Assuming an ordered state 
of the form of a coplanar spiral 



(Sj) = m[cos (p • v.i)x + sin (p • r^) 



(6) 



with p = 27r(l,0,0) or p = 27r(l, 1/2, 0), we predict by 
mean field theory that at T = a continuous second order 
transition occurs at Xc = 1/16, and the staggered magnetiza- 
tion vanishes for x > Xr like 



8^/2{x-Xc) 



(7) 



In the vicinity of the QCP, one may obtain a Landau ex- 
pansion of the effective action by standard methods ifSTI . The 
order parameters are the (real) staggered magnetizations $q 
at wavevectors 27rx, 27ry, 27rz (for a = 1, 2, 3): 



a=l,2,3 



(8) 



where the are the usual half-integer coordinates of the FCC 
sites using a unit length conventional cubic unit cell. The 
symmetry-allowed form of the effective lagrangian (in imagi- 
nary time r) is 



a 

+Sym[gi(<i>f)4 



+gi{^lf{^lf + g^^l^l^l^l\, (9) 

where "Sym" indicates symmetrization with respect to both 
wavevector (lower) and spin (upper) indices, and we have 
for simplicity neglected presumably unimportant anisotropy 
of the gradient terms. Note that the effective Hamiltonian in 
Eq. (jSjl has actually independent cubic "internal" spin symme- 
try and cubic "external" space group symmetry, which both 
constrain Eq. (|9]). 

This is an Euclidean multi-component field theory of 
standard type, which m D ^ d + \ = A space-time dimen- 
sions is in its upper critical dimension. Thus MFT is expected 
to be qualitatively correct, up to logarithmic corrections. Nu- 
merous properties of the ideal QCP follow directly. The gap 
A vanishes upon approaching the QCP from the spin orbital 
singlet phase, according to A y^a?7^^a;. Similarly, the 
Neel temperature vanishes approaching from the other side, 
Tjv ^ \/x — Xc- Other critical properties are readily obtained 
from the theory of a free relativistic scalar field, up to loga- 
rithmic corrections. A comparison with known experimental 
results is given at the end of this letter 

Having established the fundamental nature of the phase di- 
agram and QCP, we turn to a discussion of more subtle effects. 



Commensurate to incommensurate transition in the ordered 
phase — We first consider the effects of exchanges other than 
J2 in the ordered phase. Define now x = max{J, i^}/A, 
where {J, K} denotes all exchange coupling constants Ji, J2 
and Ki, K2- In the extreme limit x ^ 1, in which the spin 
orbit interaction can be neglected, one expects incommensu- 
rate magnetically and orbitally ordered ground states. If we 
assume the spins (pseudo-spins) form a coplanar spiral with 
wavevector p (q) and phase shift 9 ((f)) between the two fee 
sublattices, we obtain 8 conditions for such a configuration to 
be a classical ground state; 



|A(P)|-g: 

|A(q)H^ 
|A(p + q)| = 
|A(p-q)| = 



8K2 ■ 



= arg{A{p)); 
> = arg{A{q)); 
+ (f>^arg{A{p + q)); 
- (t> ^ arg{A{p- q)). 



(10) 



Here the complex function A(p) is defined as 

^ Px Py Pz . . Px . Py . Pz 

Alp — cos — cos — cos t sm — sm — sm — .(11) 

4 4 4 4 4 4 ^'^ 

In Eq. ( [TO] i, after trivially eliminating 0, 9, there are actu- 
ally six independent equations which can completely deter- 
mine the six real components of p and q. When some cou- 
pling constants vanish, a variety of different degenerate clas- 
sical ground states can be obtained ||20| . For example, when 
Ki — K2 — 0, Tiex reduces to the familiar J1-J2 model and 
the degenerate spin spiral momenta p's form the well-known 
spiral surface in momentum space Q El [151. 

When we turn on the spin-orbital interaction (A ^ 0), we 
expect the spin and orbital ordering to become commensurate 
with increasing x. This is because, as remarked above, TYq 
has axial cubic anisotropy, and moreover, it favors orbital or- 
der with twice the momentum of the spin spiral. A general 
spin spiral for which all spins are axially oriented satisfies 
p — ^(rii, ri2, ris), with ni,n2 and either all even inte- 
gers or all odd integers. Assuming Ji/(8J2) is not too large 
(as expected both from comparison with the structurally sim- 
ilar material MnSc2S4, and from the aforementioned neutron 
data), the states with p ~ 27r(l, 0, 0), 27r(l, ^,0) have low 
exchange energy, and are favored by the gain from the spin 
orbit interaction. Therefore we expect the ordered state to be- 
come commensurate for Xc < x < Xd (see Fig. [T]). For the 
Ji-J2-^ model, by comparing the classical energies given by 
the incommensurate spiral momenta from the spiral surface 
A(p) = and the commensurate spiral momentum, we 
find 



Xcl = J2/ACI «0.61(J2/Jl)' 



(12) 



Including non-zero Ki or K2, the critical Ad is expected to 
be somewhat smaller than the one found in Eq. ([12]). Since 
J2/J1 is expected to be large, we have Xd > I'c — 1/16. 

Excitations in the spin orbital singlet phase — For small 
X, deep in the disordered phase, one can obtain the excitation 
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FIG. 2: (color online). The low-lying spectrum along [100] direction. 
The solid curves are calculations from the small x expansion, with 
A — 25( J2 + K2)- The blue (dark) and green (light) curves have 
Ji + Ki — and Ji + Ki — J2 + K2, respectively. The red 
(dashed) curve is a schematic spectrum close to the quantum critical 
point, where it shows a Dirac-type structure at low energy. 

spectrum as an expansion in the exchange. At the leading or- 
der, we find the lowest-lying states form a triplet, with energy 

w(k) = A + (2J2 + 2X2)^cos(A-k) 

{A} 

-|(2Ji+2iri)^exp(za.k)| , (13) 

{a} 

where the a and A are summed over the 4 NN and 12 NNN 
lattice vectors, respectively. Along [100] direction, the energy 
minima of u; (which are also the global minima) are at — 
±4arccos [g^^^t^](see Fig.|2ji. The energy gap (for | Ji + 
Ki\ < 8{J2 + K2))is 

Comparison with experiments — Perhaps the strongest in- 
dication of proximity to a quantum critical point in FeSc2S4 is 
in the neutron scattering experiments by Krimmel et al 191 [TOl . 
A magnetic excitation is observed with a bandwith exceeding 
20K, in agreement with the expectations from the Curie-Weiss 
temperature Qcw — 45. IK. However, a much smaller gap 
A K, 2K is observed near the minimum, the wavevector of 
which is consistent with k — 27r(100) which we have argued 
is most natural for FeSc2S4. A gap A k, 2K was also observed 
by NMR measurements of the 1/Ti relaxation of Sc nuclei 
|8 |. Theoretically, we expect a double crossover between ac- 
tivated behavior for T ^ A to linear behavior 1 /Ti ^ T in 
the quantum critical regime A <C T ^ Tq, where Tq is a cut- 
off temperature of order Qcw^ and finally 1/Ti ^ const, for 
T ^ Tq. The low and high temperature limits are clearly ob- 
served, but the quantum critical behavior is not immediately 
apparent [8|. The uniform magnetic susceptibility remains 
large at low temperature despite the gap, which we take as 
a strong indication of the importance of spin-orbit interaction, 
i.e. Eq. (|4|i. The specific heat exhibits approximate power- 
law growth C^{T) w AT^-S for 0.2K< T <2K, with a linear 



term C„ (T) ~ jT below 0.2K and more complex behavior 
above 2K. Due to inversion disorder present in such spinels, 
we expect the very low temperature 7 term is attributable to 
two-level system defects, and the T^-^ behavior may represent 
a crossover from this to the T'^ magnetic contribution expected 
near the QCP. 

This work suggests numerous future directions for theory 
and experiment. Theoretically, the effects of fluctuations on 
the critical properties, especially with Ji ^ 0, warrant more 
detailed investigation, as do the effects of disorder - a rele- 
vant perturbation at the QCP. More theoretical studies that pre- 
dict experimental signatures are also warranted, such as possi- 
ble signs of the spin-orbital singlet in the Jahn-Teller phonon 
spectra. Experimentally, it would be most exciting to attempt 
to drive FeSc2S4 past the QCP into an ordered state. This 
could perhaps be accomplished by pressure, or with a suffi- 
ciently strong applied magnetic field. 
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